Extreme active molecular jets in L 1448C 
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(3 . ABSTRACT 

d \ The protostellar jet driven by L1448C was observed in the SiO J=8-7 and 

CO J=3-2 hnes and 350 GHz dust continuum at ~l"resolution with the Sub- 

if^ \ miUimeter Array (SMA). A narrow jet from the northern source L1448C(N) was 

observed in the SiO and the high-velocity CO. The jet consists of a chain of 

P^ I emission knots with an inter-knot spacing of ~2"(500 AU) and a semi-periodic 

velocity variation. These knots are likely to be the internal bow shocks in the 

jet beam that were formed due to the periodic variation of the ejection velocity 

with a period of ~15 — 20 yr. The innermost pair of knots, which are significant 

in the SiO map but barely seen in the CO, are located at ~1"(250 AU) from the 

j^ I central source, L1448C(N). Since the dynamical time scale for the innermost pair 

is only ~10 yr, SiO may have been formed in the protostellar wind through the 

_ gas-phase reaction, or been formed on the dust grain and directly released into 

ff^ ' the gas phase by means of shocks. It is found that the jet is extremely active with 

O' 
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a mechanical luminosity of ~7 Lq, which is comparable to the bolometric lumi- 
O ' nosity of the central source (7.5 Lq). The mass accretion rate onto the protostar 

VO ' derived from the mass- loss rate is ~10~^ Mq yr~^. Such a high mass accretion 

(^ I rate suggests that the mass and the age of the central star are 0.03-0.09 M© and 

(4-12) X 10^ yr, respectively, implying that the central star is in the very early 
stage of protostellar evolution. The low-velocity CO emission delineates two V- 
/\ • shaped shells with a common apex at L1448C(N). The kinematics of these shells 

c^ ■ are reproduced by the model of a wide opening angle wind. The co-existence 

of the highly-coUimated jets and the wide-opening angle shells can be explained 
by the "unified X-wind model" in which highly-coUimated jet components corre- 
spond to the on-axis density enhancement of the wide-opening angle wind. The 
CO J=3-2 map also revealed the second outflow driven by the southern source 
L1448C(S) located at ~8.3"(2000 AU) from L1448C(N). Although L1448C(S) is 
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brighter than L1448C(N) in the mid-IR bands, the momentum flux of the out- 
flow from L1448C(S) is two or three orders of magnitude smaller than that of 
the L1448C(N) outflow. It is hkely that the evolution of L1448C(S) has been 
strongly affected by the powerful outflow from L1448C(N). 



Subject headings: ISM: individual (L1448C) 
molecules — shock waves — stars: formation 



ISM: jets and outflows — ISM: 



1. INTRODUCTION 



It is known that several deeply-embedded young stellar objects are driving highly- 
collimated molecular outflows with an ex tremely-high ve locity (EHV) and jet-like component 
running along the axes of the lobes (e. g. . lBachilleiill996l ) . Since the EHV molecular jets hav- 
ing the terminal velocities of 50-150 km s~^ are concentrated within the narrow angles along 
the axes of the lobes and have large momenta comparable to those of the slowly moving 
(20-30 km s~^) "classical" outflows, the EHV jets are considered to be closely connected to 
the "primary jet" which is responsible for driving molecular outflows. Therefore, studying 
the physical and kinematical properties of the EHV jet will allow us to understand the prop- 
erties of the "primary jet" which provide clues to constrain the launching mechanism of the 
outflows. 



The low-luminosity (7.5 L©; iTobin et al.l (|2007| )) class source L1448C (also known as 
L1448-mm) in the Perseus molecular cloud complex (D~250 pc: e.g.. lEnoch et al.l (120061 )) is 
a spectacular example of an outflow with an EHV jet. The EHV component of this source 
was identified as the se condary peaks of the CO J=2-l spectra at ~±60 km s~^ from the 
cloud systemic velocity (JBachiller et al.lll990l ). The CO emission in the EHV range was found 
to be confined to a series of discrete clumps, called "bullets" , which are aligned along the 
axes of the outflow lobes and are symmetrically placed with respect to the central source. 
The EHV bullets are also observed in several transitions of the SiO (e.g. Bachiller et al. 



199ll : iDutrey et al.l 119971 ). Since the SiO emission has been barely d etected in quiescen t 
dark clouds because of its v ery low abundance (of the order of 10~^^; IZiurys et al.l (119891 ): 
Martin-Pintado et al.l (119921 )). the detection of SiO in EHV bullets suggests the presence of 
shocks that enhanced the SiO abundance in bullet gas by a factor of ^10^. Although the 
lower transition of SiO, i.e. J=2-l, was observed not only in the EHV bullets but also in the 
lower velocity component that delineates the tips and walls of the outflow cavities, J=5-4 
emission was confined to a pair of EHV bullets located at the closest positions to the star 
(JBachiller et al.lll99ll ). This suggests that the excitation condition of the EHV jet varies 
along the axes, and that the jet gas is highly excited in the close vicinity of the driving 
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source. Recently, higher excitation SiO up to J=ll-10 has been observed by iNisini et al. 
( 2007). Their results have revealed that the innermost pair of bullets, labeled Bl and Rl by 



Bachiller et al.l ( 119901 ). have a density of rin^ ~ 10® cm ^ and a kinetic temperature of T^in > 



500 K, which is denser and warmer than the bullets in the downstream. It is also known 
that the innermost pair of bullets, Bl and Rl are resolved into two clumps, BI-BII and 
RI-RII, respectively, in the higher resolution ( ^2") interferometric SiO J=2-l observations 
(JGuilloteau et al.l Il992l : iGirart fc AcordI l200ll ). The high resolution SiO data exhibit the 
kinematic structure of the EHV jet near the source, which shows an apparent acceleration of 
the jet up to 70 km s^^ within a r egion of 6" (^2000 AU] . The proper motion measurement 
of the SiO clumps carried out by iGirart &: AcordI (l200ll ) suggests that the outflow axis is 
inclined by 21° with respect to the plane of the sky, and therefore, the SiO clumps in the 
EHV jet are likely to be moving with absolute velocities of 180 km s~^. 

In this paper, we present the SiO J=8-7, CO J =3-2, and 350 GHz con tinuum im- 
ages obtained with the Submillimeter Arra}o (SMA) ( Ho. Moran. fc Ldl2004l ) at ~1 arc- 
second resolution, which is a factor of three higher than the previous SiO J=2-l images 
(JGuilloteau et al.lll992l : IGirart fc Acordll200ll ). A high angular resolution is crucial for study- 
ing the structure and kinematics of the jet near the base, at which the jet velocity increases 
up to the terminal velocity. In addition, higher transitions of SiO and CO in the submil- 
limeter waverange enable us to segregate the dense and warm gas in the EHV jet from the 
lower excitation gas in the cavity wall. 



2. OBSERVATIONS 

The observations of the SiO J=8-7 and CO J=3-2 lines and 350 GHz continuum 
emission were carried out with the SMA on 2006 December 5 in the extended conflguration 
and on 2006 December 25 in the compact conflguration. The two array conflgurations 
provided projected baselines ranging from 12 m to 222 m. Since the primary beam of the 
SMA antenna has a size of ~35", two pointings separated by 17"were observed in order to 
cover the EHV bullet pair closest to the central source. The receivers have two sidebands; 
the lower and upper sidebands covered the frequency ranges from 345.5 to 347.5 GHz, and 
from 355.5 to 357.5 GHz, respectively. The SiO J=8-7 and CO J=3-2 lines were observed 
simultaneously in the lower sideband. The SMA correlator divides each sideband of 2 GHz 



"'^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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bandwidth into 24 "chunks" of 104 MHz width. We used the configuration that gave 256 
channels to all chunks, which provided a uniform frequency resolution of 406.25 kHz across a 
2 GHz- wide band. The corresponding velocity resolution was 0.35 km s~^. We used Titan for 
fiux calibration, and a pair of quasars 3C84 and 3C111 for amplitude and phase calibrations. 
The flux calibration was estimated to be accurate to 25%. The band pass was calibrated by 
observing 3C273. 

The calibrated visibility data were Fourier transformed and CLEANed using the MIRIAD 
package. The velocity-channel maps of the SiO and CO were made with a velocity interval 
of 1 km s~^. The synthesized beam size of the SiO map was 0'.'96x0'.'84 with a position 
angle of —84° and that of the CO map was 0'.'96x0'.'86 with a position angle of —76° with 
uniform weighting. The rms noise level of the velocity-channel map at 1.0 km s~^ width 
was 0.15 Jy beam~^. A non-linear joint deconvolution, MO SSDI, which is ba sed on the 



CLEAN-based algorithm, and is part of the MIRIAD package (ISault et al. Ill996l ). was used 
for deconvolving the images. 

The 350 GHz continuum map was obtained by averaging the line-free chunks of both 
sidebands. To improve the signal to noise ratio, the upper and lower sidebands data 
were combined. The synthesized beam size of the map made with uniform weighting was 
0'.'93 X 0'.'83 with a position angle of —86°. The rms noise levels of the 345 GHz continuum 
maps was 6.4 mJy beani"^. 



3. RESULTS 

3.1. Continuum emission 

The 350 GHz continuum map reveals a bright compact source at the center. This source 
has a peak intensity of 352 mJy beam~^, and is surrounded by spatially extended emission. 
In addition to this bright source, there is a faint emission peak of ~5(T level at ~8'.'3 southeast 
of the center. Recent Spitzer Space Telescope observations at mid-inf rared wavelengths hay e 



resolved L1448C into two components, L1448C(N) and L1448C(S) flJorgensen et all 120061 ). 
The bright submillimeter source at the center corresponds to L1448C(N), which is considered 
to be the driving source of the highly-collimated molecular outflow and is also referred to 
as L1448-mm, and the southern faint source corresponds to L1448C (S). Dust continuum 



emission from L1448C(S) was also detected bv |j0rgensen et al.l (120071 ) at 230 GHz and 350 



GHz with the SMA, and by iMaurv et al.l (120101 1 at 107 GHz with the Plateau de Sure 
Interferometer (PdBI). 

The visibility amplitudes plot for L1448C(N) as a function of uv distance (Fig. |2]) sug- 
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gests that this source consists of two components; one is from a spatially extended envelope 
that dominates the flux at a uv distance of <50 kA, and the other is from a compact source 
that is prominent at >50 kA. The visibility amplitude profile was fit by two circular gaussian 
components; one is an extended component with a deconvolved size of ~3.6"and a flux of 
~440 mJy, and the other is a compact component with a deconvolved size of ~0.3"and a flux 
of ~330 mJy. The peak flux values of the extended and compact components correspond to 
~27 mJy and ~330 mJy beam"^, respectively, per 0'.'93 x 0'.'83 beam. As shown in Fig. la, 
the extended component appears as a bump in the northwest and a tail in the southeast (both 
of them are in the 3a level) of the central source, suggesting that this emission extends along 
the ou tflow axis. Similar emi ssion featu re along the outflow axis is also seen in the raaps o f 
3 mm JCuilloteau et al.lll992h. 2 .6 mm JSachiller et al.lllQQsh . 1.4 mm JSchoier et al.ll2004l ). 
and 1.3 mm (JMaury et al.ll201Cll ) observed with 1.5-3"resolution. Such a faint elongated fea- 
ture was clearly seen in L1157-mni , and was interpreted as the edges of the cavity excavated 
by the outflow (JGueth et al.l 119971 ). It is, therefore, possible that the extended emission in 
L1448C(N) also delineates the inner part of the envelope which is disturbed by the outflow, 
although the cavity-like structure is not clearly seen. An alternative interpretation is that 
the faint emissio n comes from an embe dded companion. A faint secondary peak seen in the 
230 GHz map of lJeirgensen et al.l (20071) implies this p ossibility. However, the position of the 
secondary peak at 230 GHz in |j0rgensen et al.l (20071) is l^^ o ffset t oward the northwest from 



that of our 350 GHz map. As discussed in |j0rgensen et al.l (120071 ) such a small difference in 
position between the 350 GHz peak and 230 GHz peak could be introduced by the extended 
emission from the envelope component sampled by different uv coverages. Since the total 
flux at 350 GHz observed wit h the SMA (extended + compact components) corresponds to 
43 % of the flux observed by iHatchell et aP J2005h at 850 /zm with SCUBA at the JCMT 
(1.737 Jy per 14"beam), it is possible that the extended component that is not sampled well 
with the SMA affects the morphology of the faint component. Although the northwestern 
bump may harbor an embedded companion, it is likely that most of the extended emission 
arises from the inner part of the envelope. Therefore, the former scenario, which attributes 
the extended component to the envelope-outflow interacting region, is more preferable in 
this case. 

The visibility amplitude profile at uv distance longer than ~50 kA implies that the com- 
pact component is not point-like but has a spatially resolved structure. In order to exclude 
the contamination from the extended envelope component, we made a map using the visibil- 
ity data with the uv distances greater than 70 kA (see Fig. lb). The synthesized beam of this 
map is 0'.'70x0'.'50 with a position angle of —87.2°. It is shown that the compact source has 
an elongation along the axis perpendicular to the outflow axis. A two dimensional Gaussian 
fit to the visibility data with uv distances larger than 70 kA yields the deconvolved major and 
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minor axes of 0.37"(90 AU) and 0.26"(65 AU), respectively, with a position angle of ~70°. 
The source position derived from the fit is a(J2000) = 3^25'^38.873^ (5(J2000)=30°44'05'.'35. 
This position agrees well (less than 0^/1) with the 3 .6 cm continuum position observed with 
a smaller beam of 0'/31x0'/27 jReipurth et alJl2002[ ). 



The 35 GHz flux of L1448C('S ) was measured to be ~60 mJy, which is c onsistent to that 
report ed bv lj0rgensen et al.l (12007 ). L1448CfS') was detected in 230 GHz by |j0rgensen et al. 



(120071) but not bylMaurv et al.l (120 lOl ) at the same frequency. The non-detection of this source 
by lMaury et al.l (120101 ) is probably because this source is located near the edge of the primary 
beam of the PdBI (~22") . If the response of their primary bea m is taken into acc o unt, t he 230 
GHz continuum source with a flux of 12.8±3 mJy detected by |j0rgensen et al.l (120071 ) could 
be below the detection limit of their observations (~8.4 mJy beam~^). Although L1448C(S) 
is bright in mid-infrared, its sub-mm flux is more than ten times weaker than L1448C(N). 
It is unlikely that such a small sub-mm flux is due to the effect of ra issing flux, because 
single-dish measurements at 450 and 350 /im (jChandler fc Riched l2000l ) showed no hint of 
the secondary component to the south of L1448C(N) even though the angular resolution 
(~8") was comparable to the separation of two sources. 



3.2. The SiO Jet 



The SiO J=8-7 emission was detected in two velocity ranges from —70 to —12 km s~^ 
(blueshifted) and from 20 to 71 km s~^ (redshifted) with respect to the systemic velocity of 
^SR~5.0 km s~^. Fig. |3] shows velocity channel maps of the SiO J=8-7 emission at 10 km 
s~^ intervals. It is shown that the SiO emission comes from the jetlike narrow region with 
its blueshifted part to the northwest and the redshifted part to the southeast of L1448C(N). 
The SiO J=8-7 jet is partially resolved along its minor axis; after deconvolution from the 
SMA beam, the width of the SiO jet is ~0.8"(~200 AU) FWHM on average. In order 
to estimate the missing flux, we have smoothed the SMA map to a resolution of 14"and 
compare d it with the SiO J =8-7 spectra observed with the James Clerk Maxwell Telescope 
(JCMT; iNisini et al. 1120071 ). It is found that 85-100 % of the single-dish flux is recovered by 
the SMA, implying that almost all the SiO J=8-7 emission arises from the narrow jet. SiO 
J=8-7 emissio n mainly comes from t he Bl and Rl "bullets" identified in the single-dish CO 
J=2-l map by lBachiller et al.l (jl990l ). The SiO jet consists of a chain of knots with a typical 
siz e scale of ~l-1.5^^ A comparison with the previous SiO J=2-l maps with ~3"resolution 



(JGuilloteau et al.lll992l : iGirart &: Acordll200ll ) reveals that the three pairs of knots close to 
L1448C(N) seen in the SiO J=8-7 map corresponds to the inner pair BI and RI in the SiO 
J=2-l map, and the two pairs in the downstream correspond to the outer pair BII and 
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RII. The innermost knot pair Bla and RIa are located within 1"(250 AU) from L1448C(N). 
The high-resolution image also shows that the SiO jet is not straight. A close-up view of 
the high velocity component (Fig. \^ shows that the jet changes its position angle from 
+15° at BII, -25° at BI, -20°at RI, to -5° at RII. The kinks between BI and BII, and RI 



and RII are also seen in th e previous SiO J=2-l maps of ~3"resolution ( JGuilloteau et al. 



19921 : iGirart fc AcordI 120011 ). However, it is more obvious in the higher resolution image. In 



addition, it is clear that the jet axes in BI and RI are also misahgned by 5°. 



3.3. CO Jet and outflow 

The CO J=3-2 emission was detected in the wide velocity range from —77 km s~^ to 
-|-79 km s^^ with respect to the systemic velocity. In the velocity ranges of AV < ±40 km 
s~^ (Fig. 4a-4d), the CO J=3-2 emission delineates two V-shaped structures open to the 
northwest and southeast with a common apex at the position of L1448C(N). The opening 
angles of the V-shape features become narrower as the velocity offset increases, suggesting 
that the CO emission in the V-shaped features comes from the limb-brightened shells. The 
largest opening angle of the shell is ~60° in the blueshifted lobe, while it is ~40° in the 
redshifted lobe. In the higher velocity ranges of AV = ±41-70 km s~^ (Fig. 4e, 4f, and 4g), 
the CO emission comes from a narrow jet-like region. The CO flux recovered by the SMA 
in each velocity range was estimate by comparing the CO spectra observed by the SMA 
with those observed by the JCMT. The SMA map was smoothed to be a 14"resolution so 
as to match with the beam of the JCMT. It is found that the recovered flux is only ~20 % 
in the lowest velocity ranges {AV = ±1-10 km s^^; Fig. 4a), and is ~50% in the next 
velocity range of AV = ±11-20 km s^^ (Fig. 4b). The edges of the shells in Fig. 4a and 
4b look very steep probably because significant amount of the CO emission from spatially 
extended component was filtered out by the interferometer. In fact, previous examples of the 
L1157 outflow and IRAS 04166+2706 outflow revealed that the edges of the bipolar cavities 
were emphasized in the maps made with the interferometer data alone, while the cavities 



were filled by the diffuse CO emis sion when the single-dish data were added (JGueth et al. 



19961 : ISantiago-Garcia et al. 1120091 ). On the other hand, in the higher velocity ranges with 
AV > ±20 km s-\ 80-100% of the CO fiux was recovered by the SMA. This suggests that 
almost all the CO J=3-2 fiux in the extremely high velocity ranges (Fig. 4e, 4f, and 4g) 
comes from the narrow jet. 

As in the case of the SiO jet, the CO J=3-2 jet also shows clumpy structure. In 
addition, most of the knots seen in the SiO J=8-7 map have their counterparts in the CO 
map. However, the innermost knots Bla and RIa, which are significant in the SiO map at 



AV = ±21-60 km s^^ are barely seen in the CO map. This is probably because most of 
the CO molecules in these knots are excited to the levels higher than J=3 because of high 
density and high temperature. Similar feature with strong SiO and weak CO in the close 
vicinity of the protostar was also ob served in the highly collimated jet in the HH211 outflow 



( jPalau et al.ll2006l : iLee et al.ll2007bl ). Another difference between the CO jet and SiO jet is 
seen between BIc and Blla, and RIc and Rlla. The CO map reveals the knot pair labeled 
BI-II and RI-II, while the SiO map shows only faint emission (Fig. [5]). In the CO J=3-2, the 
overall distribution of the EHV jet is more continuous as compared to the SiO distribution. 

As in the case of the SiO jet, the blue and red axes of the CO jet are also misaligned by 
~5°. The kinks between BI and BIT, and RI and RII are also seen in the CO jet. In addition, 
the ridge of the CO emission wiggles along the jet axis. This wiggling feature is clearly seen 
in the maps of the highest velocity ranges (Fig.E]). It is likely that each knot has been ejected 
in slightly different direction. Since the typical knot separation, 2-3"corresponds to a time 
interval of 15-20 yr, the observed jet wigghng suggests that the direction of jet ejection also 
varies in a similar time scale. On the other hand, the CO emission with lower intensity that 
is surrounding the emission ridge tends to extend linearly along the axes. The transverse 
width of the lower level CO emission component increases with a distance from the source. 
This is probably because significant part of the CO emission comes from the outflow shell 
even in the extremely high velocity ranges (see the next section). There is faint CO emission 
with the highest velocity seen ahead of the SiO jet. Since this highest velocity emission is 
spatially extended, it is likely that it arises from the highest velocity part of the shells. 



3.4. Kinematics of the jet along its axis 

Position- velocity (P-V) diagrams of the SiO and CO emission along the jet axes (the 
position angle is —25° for the blueshifted part and is —20° for the redshifted part) are shown 
in Fig. [71 Due to the change of the position angle, the outer knots Bllb and Rllb do not 
appear in these P-V diagrams. 

The velocity structure of the jet is well traced in the SiO. The jet velocity rapidly 
increases within ^l"from the star, and reaches close to its highest velocity (±65 km s~^ 
from Vgys) at ~5"from the star, which corresponds to the positions of the BIc and RIc knots. 
The velocity dispersion is extremely large (AV^~50 km s^^ at la level) at the base, while 
it narrows to AK~20 km s~^ at the positions of the BIc and RIc knots. In the positions 
further downsteam, from BIc to Blla and from RIc to Rlla, the observed radial velocity 
slightly decreases in the blueshifted side, while it increases in the redshifted side. These 
radial velocity changes are probably due to the change of inclination angle, because the 
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position angle of the jet also changes at these positions. The velocity pattern shown in 
the SiO is similar to those of the atomic jets from class I source and T Tauri stars, in 
which the low-velocity components with broad line widths are located near the base and 
the high- velocity comp onents with narrower line widths are located further from the source 



[e.g. |Pyo et al. 1120021 ). In addition to the global velocity structure, each knot shows its 
internal velocity gradient with its higher velocity in the upstream side and lower velocity in 
the downstream side. 



3.5. Kinematics of the CO outflow shells 



In the CO J=3-2, the jet shows a similar velocity pattern with a similar velocity centroid 
as the SiO. In addition to the jet, the P-V map of the CO shows extended low velocity features 
(slower than ±15 km s^^) and linear velocity features with the velocity magnitude increasing 
with the distance from the source ( "Hubble-law" ) . These Hubble- law features are obviously 
different from the jets, and are not seen in the SiO, suggesting that they are from the outflow 
shells. It should be noted that the radial velocity offsets of the linear velocity components 
become larger than 50 km s~^ at around ^10"from the star, and contaminate the CO map 
in the EHV range. 

The observed velocity pattern of the outflow shell is different from that predicted by 
the jet-driven bow shock model, in which the P-V structure shows a convex spur with hi gh- 
velocity components at the jet head (e.g. iMasson fc Chernin Ill993l : IChernin et al. Ill 9941): it 



is rather similar to the velocity pattern produced by the wide-angle wind model (IShu et al. 



199ll : IShu et al. I I2OOOI : iLi fc Shu 1119961 ). In the wide-angle wind model, an outflow shell 
consists of ambient material swept-up by a radial wind from a young star. A "Hubble-law" 
in the shell velocity is expect ed if the shell i s expanding into an ambient medium with a 
radial density profile of ocr"^ (IShu et al.lll99ll ). The shape of the shell is determined by the 
combination of the poloidal density profiles of the wind and ambient gas. It is approximately 
parabolic for an X-wind type of wide opening angle wind with an angle-dependent density 
profile of ocl/sin^^ (where is an angle measured from the axis of the flow) expanding into 
an am bient medium w ith a (xsiv?6 /r"^ density profile, which is appropriate for magnetized 
cores (JLi fc Shu Ill996l ). T herefore, we ado pted the simplified analytical model of a wind- 
driven model proposed by iLee et al.l ( I2OOOI ) to examine whether the observed morphology 
and kinematics of the CO outflow shells can be explained by means of wide-opening angle 
wind model. In the cylindrical coordinate system, the structure and velocity of the shell can 
be written as follows: 

Z = CR^, Vr = VqR, Vz = VqZ, (1) 
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where z is the distance along the outflow axis; R is the radial size of the outflow perpendicular 
to z ; C and Vq are free parameters which describe the spatial and velocity distributions of 
the outflow shell, respectively. The observed outflow shell features and the velocity patterns 
of the redshifted and blueshifted lobes were successfully reproduced by the model curves 
with C=0.8 arcsec^^and vo=5.0 km s~^ arcsec"^, and C=0.6 arcsec"^ and ?7o=5.0 km s~^ 
arcsec"^, respectively. Here the inclination angle of the outflow axis with respect to the plane 
of the sky was assumed t o be 21°, which i s deri ved from the proper motion measurement 



for the SiO J=2-l jet by iGirart &: AcordI ( l200ll ). The dynamical age of the outflow shell 



is given by I/vq and is estimated to be ~240 yr. This number is roughly consistent with 
the dynamical age of ~500 yr derived from the extent of the shell (~25"= 6300 AU), the 
mass-weighted-mean radial velocity of the shell (~22 km s~^), and the inclination angle of 
the outflow axis (~21°). The model curves that delineate the outer boundaries of the lobes 
projected onto the plane of the sky are shown in Fig. [H] on top of the contours of the outflow 
shells. The model curves of the P-V maps are shown in Fig. [71 However, this simplified wind- 
driven model has difficulty in reproducing several observed features. First, the observed CO 
intensity drops sharply at the systemic velocity and does not extend to the opposite velocity 
ranges, while the model curves on the P-V maps predict the emission at Vlsr~10 km s~^ in 
the blueshifted lobe and at ^lsr^O km s~^ in the redshifted lobe. Second, the shapes of the 
shells in the different velocities cannot be reproduced (Fig. |9]). The observed CO emission 
shows V-shaped distributions in the velocity ranges close to the systemic velocity. On the 
other hand, the model curves predict eUiptical shapes (Fig. 9a and 9b). In addition, the 
transverse width of the CO shell becomes narrower as the velocity offset increases, while 
the model predict the opposite trend. It should be noted that the CO emission near the 
protostar in Fig. 9c and 9d arises from the jet. 



3.6. Kinematics across the jet 

In order to search for the signs of jet rotation, the velocity gradient along the minor 
axes of the jet was examined. The P-V diagrams of the SiO and CO at the positions of 
innermost two pairs of knots (Fig. [TTj) show no clear velocity gradient at the positions of 
Bla, RIa, and Rib. On the other hand, there is some hint of velocity gradient in the SiO at 
the position of Bib; the southwestern side of the jet tends to be more blueshifted than the 
northeastern part. However, this velocity gradient is in the opposite sense to the rotation 
pattern of the NH.c; core with its b lueshifted part in the northeast and the redshifted part in 



the southwest (ICuriel et al.lll999l ). This suggests that the observed velocity gradient is not 



due to the rotation. 
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3.7. Physical properties of the EHV jet 

The physical parameters of the EHV jet were estimated using the CO flux measured 
in the velocity ranges of A^=50-70 km s~^, in which most of the CO flux was recovered 
by the SMA. We assumed that the CO emission in these velocity ranges is optically thin, 
and that the excitation condition of the CO follows the LTE. A fractional abundance of the 
CO and a mean atomic weight were adopted to be 10~^ and 1.41, respectively. The kinetic 
temperature of molecular gas in t he EHV jetwas derived to be 500-1 200 K from the far 



infrared lines of CO, H2O, and H2 ( Nisini et al 



1999 



the millimeter and submillimeter SiO lines (JNisini et al. 1120071 ). However, it is uncertain 



Nisini et al.ll2000l) . and >500 K from 



whether the CO (3-2) emission arises from the sar ne gas compon e nt that contributes to the 
higher transition lines in the far-infra red. In fact. iGusdorf et al.l ( l2008al ) have modeled the 
multi-transition SiO data observed by iNisini et al. I (120071 ) using their face-on C-type shock 
models, and obtained much lower temperature of 70-90 K. This implies that the most of the 
gas in the EHV jet is not so warm as >500 K. Therefore, an excitation temperature of ~100 
K was assumed for the calculation. 

The mass for each of the blueshifted and the redshifted jets is estimated to be ^10~' ^ 
M0, which is a few times higher than the bullet mass estimated by iBachiller et al.l (jl990l ) 
using an assumption of Tcx~20 K. The momentum and kinetic energy of the bipolar jet are 
estimated to be 0.3 Mq km s~^ and 5.0x10^^ erg, respectively. Here, the calculation was 
done using the de-projected jet velocity under the ass umption that the jet axis is inclined 
by 21° from the plane of the sky (JGirart &: Acordll200ll ). The dynamical timescale of the jet 
derived from the length (~ 20"= 5000 AU) and mass weighted mean (de-projected) velocity 
is estimated to be only ~150 yr. Because of the rather large mass and short timescale, the 
obtained mass loss rate is also large, ~10~^ Mq yr~^. Furthermore, the high velocity of 
~160 km s~^ (corrected for the inclination of 21°) brings extremely large momentum supply 
rate of ~2xl0~^ Mq km s~^ yr^^ and mechanical luminosity of ~26L0 for the jet. 

It should be noted that the derived total mechanical luminosity of ~26 Lq is a factor 
of 3.5 larger than the bolometric luminosity of the central source (7.5 Lq). One possible 
reason for such a discrepancy is that the mass of the jet was overestimated. Since the CO 
J=3-2 emission comes from both the jet and shell, the high velocity part of the shell may 
contaminate the CO flux that was used to calculate the mass of the jet. However, this effect 
should not be signiflcant, because the CO flux from the shell does not dominate the CO flux 
at AV >50 km s~^ velocity range. The second possible reason is the excitation temperature. 
The assumed excitation temperature of Tex~100 K is much l ower then the kinetic tempera- 



ture of >500 K derived from the far infrared measurements ( Nisini et al Ill999l : iNisini et al. 



2OOOI ) . However, the assumption with higher excitation temperature increases the mass and 
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dynamical parameters (i.e. momentum, kinetic energy, mass loss rate, momentum rate, and 
mechanical luminosity). For example, the assumption of Tex= 500 K increases the mass and 
dynamical parameters by a factor of 3.8, and brings an extremely large mechanical luminos- 
ity of ~90 Lq. This implies that the bulk of the molecular gas in the jet is not so warm 
as 500 K, and that the warm gas which contributes to the far infrared emission is not the 
major component. On the other hand, the lower excitation temperature of 40-50 K reduces 
the mass. However this effect is limited to a factor of 1.5 at the most. Another possible 
reason is the CO abundance, which was assumed to be 10~^. Th is value can be as hig h 



as 4xl0~^ in the chemical model of protostellar winds proposed by iGlassgold et al.l ( 119911 ). 
in which molecules such as CO and SiO are formed via gas-phase reactions in an initially 
atomic protostellar wind. If the higher CO abundance of 4xl0~^is adopted, the mass and 
all dynamical parameters are reduced by a factor of 4, and the mechanical luminosity of 
the jet becomes ~7 Lq, which is comparable to the bolometric luminosity of the central 
source. However, this value for the mechanical luminosity of ~7 Lq should be the lower 
limit, because it is derived under the assumption of the optically thin CO emission. The 
physical parameters derived using Tex = 100 K and CO/H2 abundance ratio of 4xl0~^ are 
given Table 1. 



3.8. Physical parameters of the outflow shell 

The mass, momentum, and kinetic energy of the outflow shell were estimated using the 
CO flux measured in the velocity ranges of ±1-40 km s~^ from the systemic velocity. The 
CO emission is assumed to be optically thin and in the LTE condition with an excitation 
temperature of 40 K, which is derived from the observed peak brightness temperature of the 
CO J=3-2 line. Since most of the gas in the shell component is considered to be the swept- 
up ambient material, the canonical CO abundance of 10~^ was used for the calculations. 
The dynamical parameters of the outflow shell are summarized in Table 2. Since signiflcant 
fraction of the CO flux is missed in the low velocity ranges. Table 2 show the parameters 
with and without correction for the effect of missing flux. The inclination angle of the 
outflow axis is assumed to be 21° as in the case of the EHV jet. The dynamical timescale 
of the outflow shell is adopted to be ~240 yr, which is derived from the modeling descried 
in the previous section. Table 2 shows that the dynamical parameters of the outflow shell 
are comparable to those of the EHV jet listed in Table 1. If the effect of the missing flux is 
corrected, the outflow shell has a momentum supply rate of ~1.8xl0~^ Mq km s~^ yr~^, and 
a mechanical luminosity of ~8.6 Lq. The mechanical lumiosity derived here is comparable 
to the bolometric luminosity of the central source and the mechanical luminosity of the 
EHV jet. It should be noted that the dynamical parameters for the redshifted outflow are 
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affected by the contamination from the L1448C(S) outflow (see the next section). However, 
the effects of L1448C(S) outflow to the dynamical parameters are not so signiflcant, because 
this outflow is seen only in the low velocity ranges. 



3.9. CO outflow from L1448C(S) 

In the CO maps in the low velocity ranges such as Fig. 4a and 4b, there is a blueshifted 
component in the redshifted lobe. This blueshifted component shows a triangle shape with 
its apex at the position of L1448C(S), and extends to the northeast direction with a posi- 
tion angle of ~40°. It is likely that this blueshifted component is related to the activity of 
L1448C(S). Fig. [11], which provides a close-up view of the L1448C(S) region, shows a red- 
shifted counterpart to the southwest of L1448C(S), a compact component at ~l"southwest 
of L1448C(S) and another extended component to the southwest of the V-shaped shell of the 
L1448C(N) outflow. This NE-SW outflow from L1448C(S) is also highly collimated. The 
opening angle of the lobes is ~40°, which is similar to the redshifted lobe of the L1448C(N) 
outflow. However, this outflow is seen only in the velocity ranges slower than 15 km s~^, 
and has no high-velocity jet-like component. There is no SiO J=8-7 emission either. 

The redshifted part of the L1448C(S) outflow overlaps with the western wall of the 
L1448C(N) outflow. At the place where the two outflows are superposed, the CO emission 
is signiflcantly enhanced and the wall of the L1448C(N) outflow lobe is bending. Therefore, 
it is possible that the two outflows are intersecting, although three dimensional geometries 
of two outflows are uncertain. 

The physical parameters of this outflow were estimated assuming that the CO emission 
is optically thin and that the excitation condition of the CO is in LTE. The fractional 
abundance of the CO was adopted to be the canonical value, 10""^, because the V-shaped 
morphology and low velocity suggest that the bulk of the L1448C(S) outflow is likely to be 
the swept-up ambient gas. A mean atomic weight of the gas was assumed to be 1.41. The 
excitation temperature was assumed to be 20 K, which is same as the rot ational temperatur e 



of the dense gas surrounding L1448C derived from the NH3 observations (jCuriel et al.lll999l ). 
The mass of the blueshifted component is estimated to be 5.3x10"'^ Mq and that of the 
redshifted component is 2.7x10"'^ Mq. Since the CO J=3-2 emission is assumed to be 
optically thin, the values derived here are the lower limits. In addition, the mass of the 
redshifted component is likely to be underestimated, because we exclude the region where 
the two outflows are superposed. It is also possible that part of the flux from the spatially 
extended component is missing. The dynamical parameters of the flow are summarized in 
Table [31 Table 3 gives two kinds of values; one is corrected for the inclination effect and the 
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other is without this correction (uncorrected). Here, the inchnation angle of the L1448C(S) 
outflow from the plane of the sky was assumed to be 32.7°, which corresponds to the mean 
inclination angle from assuming randomly oriented outflows. We also assumed that the 
outflow lobes have a size of ~20"(5000 AU). The actual size of this outflow is not clear because 
of confusion with L1448C(N) outflow and its si delobes. However, th ere is no counterpart of 
this outflow in the single-dish CO J=3-2 map ( jHatchell et al.ll2007l ). This suggests that the 
spatial extent of this outflow is not large, or the spatially extended component is fainter than 
the sensitivity limit of the single-dish measurement. In the later case, the faint component 
would not make significant contribution to the dynamical parameters even though its spatial 
extension is large. 



4. DISCUSSION 



4.1. Compact disk around L1448C(N) 



The compact component of dust continuum emission is partially resolved by the 0'.'7x0'.'5 
beam. The observed structure is elongated perpendicular to the outflow axis, suggesting 
that this component traces the disk surrounding the protostar. If the beam-deconvolved 
major axis, ~90 AU, represents the diameter of the disk, the disk size of L1448C(' N) is 
compara ble to those around the youngest pro t ostars such as H H211 (JLee et al.ll2007bf ) and 
HH212 JCodella et al.l l2007t llee et al I [2007al : llee et al.ll2008l ). If the measured flux (330 
mJy) comes from the region with 0'.'37x0'.'26 size, the brightness temperature corresponds 
to ~35 K. The spectral energy distribution (SED) of the compact source from 8.3 GHz to 
350 GHz is shown in Fig. [121 The measured flux densities are fit by a single power law 
with a spectra l index a of 1 . 98 (so lid line in Fig [T2|) . which agrees with the previous result 
[a = 1.84) of ISchoier et al.l ( 120041 ). This spectral index is smaller than a = 3.4 derived 
from the photometric broadband mea surement including the contribution of the larger scale 
envelope component (JFroebrichI |2005[ ) . Since the emission at cm wavelengths is likely to be 
the free-free emission from shock-ionized gas, the contribution of the free-free c omponent at 
mm a nd sub-mm wave ranges was estimated by usii ig the equations gi ven by ICuriel et al. 
( Il990l ) . Using the parameters of the stellar wind in ICuriel et al.l ( ll990l ) , the flux densities 
of the free-free emission at mm and sub-mm wave ranges were estimated to be less than 1 
mJy (dotted line in Fig [T2|) . Therefore, the small index number is unlikely to be due to 
the contribution of the free-free component. The power law fit assumes that the emission is 
optically thin and that the Rayleigh- Jeans approximation is valid. However, the observed 
index a ~2 is close to the value for the blackbody radiation, suggesting that the emission is 
optically thick. In addition, a brightness temperature of 35 K suggests that the Rayleigh- 
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Jeans approximation is not applicable in the mm and sub-mm wave ranges. Therefore, we 
applied an optically thick fit without Ray leigh- Jeans approximation using the formula, 



Su = QsB,y[l - exp(-rj,)]. 



(2) 



Where Si, is the fiux density, Qg is the source size, Bj, is the Planck function, and t^, is the 
dust optical depth, which is assumed to follow the power low, Ty(xv^ . Assuming that the 
source size is 0'.'37x0'.'26 and a dust temperature of 40 K, the SED fit provides /3 = 1.3 and 
T350GHZ = 7.5 (dash-dotted line in Fig [12]). The average optical depth at 350 GHz for a disk 
of mass (gas+dust) M^ and radius Rd is given by 



< T350GHZ >= 



0.5 



M, 



cos^y \q.iMqJ ViooAuy 



(3) 



where 9 is the disk inclination angle to the line of sight (jjgirgensen et al.l 120071 ). Using this 
relation, the mass of the disk with tssoghz = 7.5, Rd = 45 AU, and 6 = 69° (assuming that 
the disk is perpendicular to the jet axis) is estimated to be 0.11 Mq. The mass derived here 
is approximately twice as large as the lower limit of 0.047 Mq, which is derived under the 
assumption of optically thin emission. 



4.2. Stellar mass loss rate and its implication to the protostellar evolution 

Since protostellar jet is considered to be closely linked to the mass accretion, the stellar 
mass-loss rate gives us a rough estimate of the mass accretion rate onto the star. Theoretical 
estimate for the ratio of mass outfiow to mass accretion rate {Mout/Macc) is ~l/3 for an X- 
wind type magneto-centrifugal wind (e.g. IShu et all 1 19941 ). If the Mout/M^cc ratio is assumed 
to be ~0.3, the total mass-loss rate (blue -|- red) derived from the CO fiux, 2.4x10"^ Mq 
yr~^, gives us the mass accretion rate of 8x 10~^ Mq yr^^. In spite of rather high accretion 
rate, the observed bolometric luminosity is only 7.5 Lq, suggesting that the mass of the 
central star is still very small. If most of the observed bolometric luminosity is released by 
means of accretion, the mass of the central star is calculated by using the relation, M^, = 
Lacc-R*/GMacc, wherc M* is the mass of the central star. Lace is the accretion luminosity, 
and Mace is the mass accretion rate onto the protostar. The radius of the protostar, i?*, 
is considered to be ~1 -R^r ^ in the earlies t evolutionary stage with very low mass and ~3 
Rq in the later stage (e.g. IStahler 1119881 ). The mass of the central star is estimated to be 
0.03-0.09 Mq for an accretion luminosity of 7.5 Lq. With a stellar mass of 0.03-0.09 Mq, 
the Keplerian velocity at the surface of the protostar becomes ~80 km s^^. In this case, 
the jet velocity to Keplerian velocity r atio becomes ~2, which is reasona ble if the jet is 
launched by magneto-centrifugal force (IShu et all 1 19941 : iPudriz et al. II2007I ). If a constant 
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mass accretion rate of 8x 10^^ Mq yr~^ is assumed, the age of the central star is estimated to 
be (4-12) X 10^ yr. The timesc ale derived here is con sistent to the kinematic age of the larger 
scale outflow of ~0.3 pc scale ( iBachiller et al.lll990l ). However, morphology of the EHV jets 
implies that the mass accretion was variable. The highly-colimated EHV jets terminate at 
~20"from the source, suggesting that L1448C(N) experienced lower activity phase in the 
past and enhanced its activity significantly in the last ~150 yr. 



4.3. Clumpy structure in the L1448C(N) jet 

High resolution S iO and CO maps show that the EHV bullets Bl and Rl identified by 
Bachiller et al.l ( 1l990l ) consist of chains of knots. If the BI-II and RI-II knots are included, 
the knots are aligned with almost equal intervals of ~2"(500 AU). Simil ar knotty structure 



with semi-regular intervals is also se en in the SiO and CO jets in HH211 (iHirano et al. 



Palau et al. 



20061 : iLeeet al.ll2007bh and HH212 fICodella et al.l 120071 : iLee et alll2007al ). In 



2006 



the case of HH211 and HH212, the knots seen in the SiO and CO have their counterparts in 
the near infrared H2 emission except the innermost knots pairs that were highly obscured. 
In the case of the L1448C(N) outflow, H2 emission kno ts are seen only i n the norther n 
blueshifted side and not in the southern redshifted side ( jDavis et al.lll994J : lEisloffelll2000l ). 
This is probably because the axis of the L1448C(N) jet is inclined from the plane of the sky 
and the near infrared emissio n in the southern part is obscured by the dense gas envelope 
traced by the NH3 emission ( ICuriel et al.lll990l ). In the norther n side, the morphology of 
the SiO jet coincides well with that of the H2 jet ( lEisloffelll2000l ). which also shows a kink 
between BI and BII components. Therefore it is likely that t he knots in the L1448 jet are the 



internal bow shocks in the jet beam as in the cases of HH211(lHirano et al.ll2006l : iPalau et al. 



20061 : iLee et al.ll2007bl ) and HH212 fICodella et al.l 120071 : iLee et alll2007al ). In fact, some of 
the SiO knots are partially resolved in the transverse direction, and the Rll-a knot shows 
an arc-shaped structure typical of a bow shock (Fig. 3f). The SiO emission is weak at the 
positions of the BI-II and RI-II knots, suggesting that the shocks at these positions are rather 
weak as compared to the other knot positions. Since the jet is deflected at the positions of 
BI-II and RI-II knots, it is likely that the jet material there is impacting less dense material 
surrounding the jet beam. 



4.4. Jet bending 

As shown in Fig. 5, the blue part and the red part of the jet are misaligned by ~5° 
and forming a C-shaped structure bending toward the west. Such a C-shaped bending 
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of th e jet could be due to th e Lorenz forces between the jet and interstellar magnetic 



field (Fendt &: Zinneckei 



( JFendt fc Zinneckei 



1998 



I 



9981) . the orbital motion of the jet source in a binary system 



Masciadri &: Raga 1 12002| ) , or dynamical pressure from external 



medium (JFendt &: Zinneckerlll998l ). In the case of the Lorenz forces, a C-shaped bending 
i s expected if the poloid al current in the jet and counter jet flows in the same direction 
( jFendt fc Zinneckerlll998l ). However, this mechanism is difficult to account for the observed 
bending of the L1448C(N) jet, because typical interstellar magnetic field with several tens 
of microgauss is not strong enough to bend the jet beam with a density of >10^ cm~^. 

If the C-shaped bending is produced by the orbital motion of a binary s ystem, the orbital 
radius and orbital velocity can be estimated by using the analytical model of lMasciadri fc Raga 
( |2002| ). Here the jet is assumed to be ejected at a velocity of Vj from one of the binary proto- 
stars in a circular orbit of radius tq and orbital velocity vq. The ^:- axis is parallel to the orbita l 
rotation axis, and z=0 is the orbital plane. As shown in Fig. 3 of lMasciadri fc Raga I ( 12002| ). 
the defiection angle a of the jet beam near the source is approximated by the x=K,z/cosi 
line, where k=vq/vj and i is the inclination angle of the orbital axis with respect to the plane 
of the sky. In the case of the L1448C(N) jet, the defiection angle a is estimated to be 2.5°, 
which corresponds to the half of the misalignment angle. Using the jet velocity of ~160 km 
s^^, the orbital velocity is calculated to be 6.5 km s~^. Since the mass of the protostar with 
jet is only 0.03-0.09 M©, the total mass of the binary system is considered to be less than 
0.18 Mq. Therefore, the radius and period of the orbital morion are estimated to be smaller 
than 4.2 AU and 20 yr, respectively. However, such a short period orbital motion cannot 
account for the observed C-shaped pattern, because the C-shaped bending is seen in the BI 
and RI parts of the jet with a length of ~2000 AU with a dynamical time scale of 47 yr. In 
order to produce the C-shaped pattern with the orbital motion, the orbital period needs to 
be longer than twice of the dynamical time scale. 

In the case of dynamical pressure of external medium, ambient gas with n(H2)<10^ 
cm~^ cannot account for bending the jet with a density of >10^ cm~^, unless the protostar 
is moving with a velocity that is comparable to the jet velocity. On the other hand, the 
dynamical pressure caused by the outfio w from the nea r by pr otostar, L1448N, cannot be 
ruled out. As shown in the CO map of iBachiller et al.l (Il990l ). the redshifted lobe of the 
L1448N outfiow overlaps with the blue lobe of the L1448C(N) outfiow. The interaction 



betwe en the two outfiows from L1448C(N) and L1448N has been suggested by IBachiller et al. 



(119951 ). because the large scale outfiow from L1448C(N) shows a considerable bending at the 
place where the two outfiows are overlapping. Sin ce the redshifted emi ssion from L1448N 
outfiow reaches close to the position of L1448C(N) (IBachiller et al.lll995l ). it is possible that 
the jet from L1448C(N) is propagating under the infiuence of the L1448N outfiow. In this 
case, the dynamical pressure from L1448N outfiow acts from north to south, and defiects the 
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jet beams to the west if they were ejected to the the northwest and southeast directions. 



4.5. Deflection and wiggling of the jet 

In addition to the C-shaped bending, the jet is also deflected toward the east by ~40° 
at the position of Bl-11 and toward the south by ~15° at the Rl-Il position. Since both 
sides of the jet are deflected at almost same distance from the central star, the jet deflection 
is likely to be caused by some variability intrinsic to the driving source rather than by 
external perturbation. The observed morphology is similar to the S-shaped point-reflection 
symmetric pattern that is expected if the disk is precessing or wobbling. Although the jet 
is not exactly the S-shape but asymmetric in deflection angle, this is probably because of 
the projection effect. In a binary protostellar system with a disk misaligned with the orbital 
plane of the binary, the disk wobbles with a period appro ximately half of t he binary orbital 



period and precesses with a period of ~20 orbital period ( iBate et al.l 120001 ). S-shaped point 
symmetry will be observed if the precession or wobbling time scale is longer than four times 
of the dynamical timescale. Since the jet deflection occurs at BI-II and RI-II the time scale 
of which is ~50 yr, the time scale of the precession or wobbling should be longer than ~200 
yr. Therefore, if the deflection is due to the precession, the lower limit of the binary orbital 
period is ~10 yr. Since this orbital period of ~10 yr is comparable to the period of small 
scale wiggling shown in Fig. 6, 15-20 yr, this orbital motion can also explain the wiggling 
feature. If the binary system consists of equal mass protostars with 0.03-0.09 M©, the orbital 
radius is estimated to be 2.4-4.2 AU. On the other hand, if the jet deflection is due to the 
wobbling of the disk, the orbital period and the separation of the binary are estimated to 
be ~400 yr and 30 AU, respectively. Since the estimated separation of the binary is smaller 
than the size of the disk observed with the 350 GHz continuum emission, it is possible that 
the observed 90 AU scale disk harbors two sources separated by 60 AU. However, the binary 
with a separation of >60 AU cannot account for the small scale wiggling feature. 



4.6. Velocity variation of the jet 

The P-V diagram of the SiO shows that the velocity of the jet varies semi-periodically. 
The velocity variation is more obvious in Fig. [131 which plots the velocity centroid of the 
SiO emission in the redshifted part of the jet as a function of the distance from L1448C(N). 
It is shown that the typical amplitude of the variation in velocity centroid is ~7 km s~^. 
Such a velocity variation is expected if the jet is precessing, the jet is launc hed from an 



orbiting object, or the ejection velocity itself varies as a function of time (e.g. ISmith et al. 
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19971 ). The period of the velocity variation estimated from the de-projected jet velocity and 
knot separation is ~15-20 yr. Since this time scale is much shorter than the precession time 
scale that is estimated to be 200 yr, it is unlikely that the velocity variation is caused by 
the precession of the jet. On the other hand, the orbital motion of the driving source can 
account for the velocity variation; the binary system with an orbital period of ~15-20 yr, an 
orbital radius of 2.4-4.2 AU, a total mass of 0.06-0.18 Mq has an orbital velocity of ~4.7-6.2 
km s~^, which is comparable to the amplitude of the radial velocity variation. However, the 
orbital motion cannot explain the relation between the SiO intensity and velocity gradient. 
As shown in the P-V map (Fig. [7]), each SiO knot has its higher velocity in the upstream 
side and lower velocity in the downstream side. The opposite velocity gradient is always 
seen in the faint emission between the knots. Such a structure is more likely to be formed 
by means of periodic variation of the ejection velocity. In such a case, the SiO knots are 
considered to be forme d as the fast moving material plunge into t he slow moving material 
in the downstream (e.g. IStone fc Normanlll993l : ISuttner et al.lll997l ). The periodic variation 
in the ejection velocity is probably due to the modulation of mass accretion by means of 
companion. In such a case, the variation amplitude of the jet velocity corrected for the 
inclination is calculated to be ~20 km s^^. This velocity amplitude corresponds to the 
shock velocity, which is consistent to the velocity of C-typ e shocks that can account for 



the e xcitation conditions of the far infrared molecular lines (JNisini et al Ill999l : iNisini et al. 



20001 ) . Similar velocity gradients in the knots with the faster part in the upstr eam side and 



the s lower part in the downstream was also observed in the o ptical jet of HHlll (Raga et al 



20021) and in the CO and SiO jets from IRAS 04166+2706 flSantiago-Garcia et al. 1120091 ) 



4.7. Driving mechanism of the CO outflow 



The P-V diagram of the CO J=3-2 along the axis (Fig. [7]) exhibits two kinematic 
components, i.e. the EHV jet with an almost constant velocity and the outflow shell with 
a parabolic velocity pattern. Although highly-collimated jet is clearly seen in both SiO and 
CO, the parabolic velocity pattern seen in the outflow shells is reproduced by the wind- 
driven model. Therefore, the observational results require a wide-opening angle wind and a 
coUimated jet at the same time. 

One possible m echanism to e xplain the observed jet+shell structure is the "unifled 
model" proposed by IShang et al.l (120061 ) . in which highly-collimated jet component is ex- 
plained as an on-axis density enhancement of the X-wind type of wide opening angle wind 
launched magnetocentrifugally. In this model, the jet along the axis corresponds to the 
densest part of the primary wind, and the shell is mostly consisted of the swept-up ambient 



-20- 



material. Their numerical model successfully reproduced the structure of a dense and narrow 
jet surrounded by a conic al shell. The o t her n aodel s that can explain the two components 
structure are proposed bvlMachida et al.l (120081 ) and lBanerjee fc Pudritzl (J2006|). The model 
proposed by iMachida et al.l (120081 ) predicts two distinct flows from the adiabatic core and 
the protostar. The flow from the adiabatic core driven by the magnetocentrifugal mecha- 
nism has a low velocity and a wide opening angle, while the flow from the protostar, which 
is mainly driven by the magnetic pressure gradient fo rce, has a high velo c ity an d is well 
collimated. On the other hand, a model proposed by iBanerjee fc Pudritzl ( l2006l ) predicts 
the structure with the jet powered by magnetocentrifugal force enclosed by the large-scale 
outflow driven by the magnetic pressure. Although these two models reproduce the jet+shell 
structure similar to the observational results, the vel ocities of the s h ells p redicted in these 
model are rather small (~5 km s~^ for the model of iMachida et al.l (120081 )) because of the 
shallow gravitational potential at the launching point. In the case of the L1448C(N) outflow, 
the terminal velocity of the outflow shellreaches to AV^~±70 km s~^ without inclination cor- 
rection, which is comparable to the velocity of the EHV jet. In order to launch such a high 
velocity wind, the launching point of this wind should be close to the launching point of the 
EHV jet. Therefore, the models with two components launched from two different regions do 
not explain the jet-l-sh ell structure in the L1448C(N) outflow. In the case of the X-wind with 
density stratification (IShang et al.l 120061 ) . the rather high velocity in the shell component is 
naturally explained because the shell is driven by the high-velocity primary wind launched 
at the same region as the EHV jet. 



4.8. Origin of the SiO in the jet 



It is considered that the SiO molecules observed in jets and outflows are formed as 
a consequence of grain sputtering in a C-shock rel easing Si-bearing rnaterial into the gas 
phase, followed by t he reaction with O and OH (ISchilke et al.l Il997l : ICaselli et al.l Il997l : 
Gusdorf et al. l2008ar ). The multi-transition SiO lines from the L1448C(N) jet observed by 



Nisini et al. 



(120071) h ave been successfully modeled by the steady-state C-shock model of 
Gusdorf et al.l (l2008al ) with a pre-shock density of ~10^ cm~^ and a shock velocity of 30- 



45 km s~^. However, the conversion of Si into SiO is initially rather slow in their models, 
and the predicted SiO line emission predominantly arises from postshock gas >100 yr after 
the passage of shocks. Since most of the SiO knots in the L1448C(N) jet have dynamical 
tir nescale shorter than t he SiO formation timescale, the steady-state C-type shock models 
of iGusdorf et al.l ( l2008al ) does not account for the SiO in the knots close to the central star, 
especially in the innermost knot pair with extremely short time scale less than 10 yr. One 
possible explanation is that the SiO molecules existed on the grain mantles and are released 
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into the gas-phase by means of shocks as suggested by iGusdorf et ahl (l2008bl ). 



A nother possibihty is the form ation of SiO in high density primary jet (IShang et al. 



20061 ). iGlassgold et ahl (119891 Il99l[ ) studied the formation of molecules in protostellar winds, 
which are originally neutral atomic, and found that significant quantities of SiO can be 
formed quickly in the close vicinity (<0.1 AU) of the central star if the mass- loss rate is high 
(>10~^ Mq). Since the mass-loss rate of the L1448C(N) jet is high enough, this scenario 
of in situ formation also ca n be the origin of the SiO in the jet. The chemical models of 
Glassgold et al.l ( 1l989l . Il99ll ) also predict significant amount of CO synthesized in the winds; 
the CO abundance reaches an equilibrium value of 4xl0~^ under the conditions in which 
observable amount of SiO is formed. The morphological and kinematical similarity of the 
CO and SiO jets supports the idea that both CO and SiO are formed in the protostellar 
wind. 



4.9. Properties of L1448C(S) 



The secondary source L1448C(S) is located at ~8'.'3 (2000 AU) southeast of L1448C(N). 
The 350 GHz continuum flux of this source is ~60 mJy, which is five times smaller than 
the flux from L1448C(N). If an optically thin condition is assumed, the mass of the cir- 
cumstellar material surrounding L1 448C(S) is estimated from t he observed flux using the 
dust mass opaci t y of 1 .75 cm^ g~^ (jOssenkopf fc Henninglll994l ) and the formula given by 
J0rgensen et al.l (120071 ). For a dust temperature of ~40 K, the estimated mass in the opti- 
cally thin limit is 8.6x10^^ Mq. Since L1448C(S) is associated with a molecular outflow, 
it is highly probable that this source is a pro tostar rather than a mere dust c l ump a t the 
cavity wall as claimed by iMaury et al.l (120101 ). The NH3 data of ICuriel et al.l (119991 ) sug- 
gest that L1448C(S) is formed in the same dense core as L1448C(N). However, the SED 
of L1448C(S) is significantly different from that of L1448C(N). In table 4, the broad-band 
spectra of L1448C(N) and L1448C(S) measured at different wavelengths are listed. In the 
near-infrared, L1448C (S) is much d imme r than L1448C(N); only L1448C(N) appeared in 
the Kg band image of iTobin et al.l (120071 ). On the contrary, in the mid-inrfrared at three 
IRAC bands, band 2 (4.5 /im), band 3 (5.8 /im), and band 4 (8.0 /im), L1448C(S) becomes 
brighter than L1448C(N); especially in the bands 3 and 4, L1448C(S) is more than six 
times brighter tha r i L1448C(N). In t he MIPS 24 /im image, L1448C(S) is also seen clearly 
jTobinet al.ll2007l : iRebuU et alll2007h . The flux from L1448C(S) at 24 /xm looks similar to 
that from L1448C(N), although the accurate flux value of each source is not easy to measure 
because of the confusion. In the sub-mm and mm wavebands, L1448C(S) is much weaker 
than L1448C(N). The masses of the circumstellar material surrounding L1448C(S), ~0.01 
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Mq, is approximately 10 times smaller than that of L1448C(N), ~0.1 Mq. Due to the 
small amount of circumstellar material, the central star of L1448C(S) is likely to be less 
obscured in the mid-infrared as compared to that of L1448C(N) enshrouded by the thick 
cocoon. The outflow activities in two sources are also different significantly. The CO out- 
flow from L1448C(S) is compact and substantially weaker than the L1448C(N) outflow. The 
momentum flux of the L1448C(S) outflow is only ~10~^ Mq km s~^ yr^^, which is two or 
three orders of magnitude smaller t han that of the L 1 448C (N) outfow and is comparable to 
those of class I outflows studied by iBontemps et al.l ( 1l996l ). In addition, there is no EHV 
component nor SiO emission associated with the outflow from L1448C(S). 

The small amount of circumstellar material of less than 0.01 Mq suggests that L1448C(S) 
may have accumulated most of its circumstellar mass. The less energetic outflow also support 
the idea that L1448C(S) has a nature close to class I rather than class 0. These results imply 
that two sources, L1448C(N) and L1448C(S) are formed in different epochs in the same dense 
core. Another possibility is the effect of the L1448C(N) outflow. Since L1448C(S) is located 
at the same line of sight as the L1448C(N) outflow, it is possible that the outflowing gas 
has stripped away the dense gas surrounding L1448C(S). Although the three-dimensional 
geometries of the sources and outflows a re not clear, the hi gh rotational temperature of NH3 
observed at the position of L1448C(S) ( ICuriel et al.lll999l ) suggests the possibility that the 
gas around L1448C(S) is heated by the interaction with the jet from L1448C(N). In such 
a case, the apparent age of L1448C(S) is older than that of L1448C(N) simply because the 
amount of material left around L1448C(S) is smaller than that around L1448C(N). The effect 
of outflow is also proposed to explain the difference of the apaprent evolutionary stage of 
protostellar pair L144 8N('A) and L1448N(B), which is located at ~75"northwest of L1448C 
Jo 'Linger etaPbooel ). 



The t hird s cenario is the disintegration of an unstable multiple system as proposed by 



ReipurthI (I2OOOI ). Since nonhierarchical triple systems are unstable, they break up, ejecting 
the lightest member, while the remaining two members form a close binary system with 
a highly eccentric orbit. In this scenario, the disks around escaping stars will be highly 
truncated; the typical disk size is expected to be around half of the distance between the 
stars in the close triple encounter. If L1448C(S) is the escaping member, the small amount 
of its circumstellar material can be explained by means of disk truncation. This scenario 
also implies that L1448C(N) is a close binary system. The observed deflection, wiggling, 
and periodic velocity variation of the jet suggest the possibility that L1448C(N) is a close 
binary system with an orbital radius of ~2-4 AU. Therefore, it is possible that such a close 
binary system was formed by means of disintegration of a triple system. In order to assess 
this scenario, kinematic iii f ormat ion of L1448C(S) becomes important. Although previous 
NH3 results of ICuriel et al.l ( 119991 ) did not show peculiar motions in the dense core containing 
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L1448C(N) and L1448C(S), a detailed study with higher angular resolution would be helpful. 



5. CONCLUSIONS 

The central region of the highly-colimated molecular outflow driven by L1448C was 
mapped in the SiO J=8-7, CO J=3-2, and 350 GHz continuum emission with the SMA at 
~1 arcsecond resolution. Our main conclusions are the following: 

1. The 350 GHz continuum emission was detected from two Spitzer sources L1448C(N) 
and L1448C(S). The continuum emission from L1448C(N) consists of an extended com- 
ponent and a compact component. The compact component is elongated perpendicular 
to the outflow axis, and is likely to be a circumstellar disk with a size of ~90 AU. The 
spectral index of this compact component derived from the data from 86 GHz to 350 
GHz is a~2, suggesting the possibility that the continuum emission is optically thick 
at 350 GHz. The mass of the disk is estimated to be ~0.11 Mq. 

2. The continuum flux from L1448C(S) is ~60 mJy, which is ~10 times lower than the 
flux from L1448C(N), although L1448C(S) is brighter than L1448C(N) in the mid 
infrared wavebands. The mass of the circumstellar material surrounding L1448C(S) is 



estimated to be 8.6x10 ^ Mp 



3. A narrow jet from L1448C(N) along the outflow axis was observed in the SiO and the 
high- velocity CO. The width of the jet measured in the SiO images is ~200 AU FWHM 
on average. The jet consists of a chain of emission knots with an inter-knot spacing of 
~500 AU. It is likely that the knots in the L1448 jet are the internal bow shocks in 
the jet beam. 

4. The dynamical timescale of the innermost pair of knots, which are significant in the 
SiO but barely seen in the CO, is only ~10 yr. It is likely that the SiO may have 
been formed quickly in the protostellar wind through the gas-phase reaction, or been 
formed on the dust grain and directly released into the gas phase by means of shocks. 

5. The low velocity CO emission delineates two V-shaped shells with a common apex 
at L1448C(N). The kinematics of this shell component is reproduced by the model 
of wide opening angle wind. Therefore, the outflow from L1448C(N) consists of both 
highly-collimated jets and shells produced by wide-opening angle wind. The observed 



jet+shell structure can be explained by the "unified model" proposed by IShang et al. 



( 120061 ). in which highly-collimated jet components are explained as an on-axis density 



enhancement of the X-wind type of wide opening angle wind. 
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6. The Jet from L1448C(N) is extremely active with a momentum supply rate of ~5xl0~^ 
Mq km s~^ yr~^ and a mechanical luminosity of ~7 Lq. The mass accretion rate 
derived from the mass loss rate is ~10~^ Mq yr~^. Such a high mass-accretion rate 
and a rather low bolo metric luminosity of the central source, 7.5 Lq, imply that the 
central protostar is still in the very early phase of its evolution with a mass of 0.03-0.09 
Mq and a dynamical age of (4-12) x 10^ yr. 

7. The blue part and the red part of the jet are misaligned by ~5° and forming a C-shaped 
bending toward the west. The possible origin of this bending is the dynamical pressure 
caused by the outflow from the nearby protostar, L1448N. 

8. The jet is deflected toward the east in the blueshifted part and toward the south in 
the redshifted part. In addition, the jet is wiggling with a period of ~15-20 yr. The 
deflection and wiggling of the jet can be explained if the driving source is a member of 
the binary system with an orbital radius of 2-4 AU. 

9. The jet shows a semi-periodic variation in radial velocity with an amplitude of ~7 km 
s~^. Each SiO knot has its higher velocity in the upstream side and lower velocity 
in the downstream side. The opposite velocity gradient is seen in the faint emission 
between the knots. It is likely that the ejection velocity varies periodically by means 
of modulation of mass accretion. 

10. The bipolar outflow in the NE-SW direction centered at L1448C(S) was discovered in 
the CO J=3-2. This provides a strong evidence that L1448C(S) is a protostar. The 
momentum flux of this outflow is only ~10~^ Mq km s^^ yr^^, which is two to three 
orders of magnitude smaller than that of L1448C(N) outflow, and is comparable to 
those of class I outflow. 

11. L1448C(S) is surrounded by a rather small amount of circumstellar material of less 
than 0.01 Mq and is powering a less energetic outflow, suggesting that this source has 
a nature close to class I rather than class 0, even though this source is formed in the 
same dense core as L1448C(N). One possible scenario to explain this dichotomy is the 
effect of the outflow from L1448C(N); signiflcant amount of material in the envelope 
surrounding L1448C(S) might have been stripped off by the powerful outflow from 
L1448C(N). Another possibility is the disintegration of an unstable multiple system, 
in which L1448C(S) is an escaping member with a truncated disk. 

We wish to thank all the SMA staff in Hawaii, Cambridge, and Taipei for their enthu- 
siastic help during these observations. N. Hirano thanks M. Machida and S. Inutsuka for 
fruitful discussion. N. Hirano is supported by NSC grant 96-2112-M-001-023. 
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Table 1. Dynamical parameters of the EHV jet from L1448C(N) 



Parameters 


Blue 


Red 


Mass (Mq) 


2.3x10-^ 


2.2x10"^ 


Momentum (Mq km s^"*^) 


0.038 


0.037 


Kinetic energy (erg) 


6.2x10^3 


6.2x10^3 


Mean velocity (km s^-^) 


160.4 


164.6 


Dynamical timescale (yr) 


160 


150 


Mass loss rate"" {Mq yr"^) 


1.2x10-^ 


1.2x10-^ 


Momentum supply rate (Mq km s^^ yr^^) 


2.3x10^^ 


2.5 xlO"^ 


Mechanical luminosity (Lq) 


3.2 


3.4 



^The primary jet velocity is assumed to be 200 km s 



-30- 



Table 2. Dynamical parameters of the L1448C(N) outflow shell 



Blue Red 



Parameters uncorrected^ corrected^ uncorrected*^ corrected^ 

Mass (M©) 1.9x10--'^ S.SxlQ-^ S.OxlQ-^ 8.6x10-^ 

Momentum (M© km s~^) 0.10 0.18 0.15 0.26 

Kinetic energy (erg) 8.1x10^^ i.ixlO^^ i.OxlO^^ 1.4x10^^ 

Momentum supply rate^ (Mq km s"^ yr^i) 4.1x10"^ 7.4x10"^ 6.1x10"^ 10.9 xlO~^ 

Mechanical luminosity'^ (Lq) 2.8 3.8 3.5 4.8 

*The effect of missing flux is not corrected. 

^The effect of missing flux is corrected. 

'^The dynamical time scale is assumed to be 240 yr. 
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Table 3. Dynamical parameters of the L1448C(S) outflow 



Blue lobe Red lobe 



Parameters i=32.7° uncorrected 2=32.7° uncorrected 

Mass (Mq) 5.3x10-^ 2.7x10-^ 

Momentum (M0 km s"^) 5.9x10^3 3.2xl0--'^ 2.0 xlO'^ l.lxlQ-^ 

Kinetic energy (erg) 8.7x10^^ 2.6x10^^ 1.7x10^^ 5.0x10^° 

Dynamical timescale (yr) 2.5 xlO^ 3.9 xlO^ 3.9 xlO^ 6.1x10^ 

Mass loss rate'' (M© yr-i) 2.4x10"^ 8.2x10-^ S.OxlQ-^ 1.7x10-^ 

Momentum supply rate (Mq km s^^r-i) 2.4x10"^ 8.2x10"^ 5.0x10"'^ 1.7x10"'^ 

Mechanical luminosity (Lq) 2.9x10"^ 5.5x10"^ 3.6x10"^ 6.7x10"^ 

''The primary wind speed is assumed to be 100 km s~^. Then we derived the mass loss rate by 
assuming momentum conservation of the the entrained CO gas and the primary wind. 
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Table 4. Broad-band spectra of L1448C(N) and L1448C(S) 





Flux (i 


iiJy) 




Wavelength 


L1448C(N) 


L1448C(S) 


Ref. 


1.65 /xm 


0.144±0.016 


— 


1 


2.1 /im 


1.07±0.13 


— 


1 


3.6 yum 


4.2±0.4 


2.5±0.08 


2 


4.5 /xm 


13.5±0.8 


23.6±0.6 


2 


5.8 //m 


11.7±0.8 


77.2±0.8 


2 


8.0 /im 


19.6±1.2 


123±1 


2 


860 /im 


330±6 


60±6 


3 


870 /im 


370±8 


44±8 


4 


1300 /im 


120±3 


13±3 


4 



References. — (1) Tobin et al. 2007; (2) 
J0rgensen et al. 2006; (3) This work; (4) 
J0rgensen et al. 2007 
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Fig. 1. — 350 GHz continuum emission from the central part of L1448-C. a) Map made using 
all the available visibility data. The contours start at 3cr and are drawn at 6, 12, 18, 27, 36, 
45, and 54cr. The la level is 6.42 mJy beam~^. b) Map made using the visibility data with 
the uv distance greater than 70 kA. The contours start at 3a and are drawn at 6, 12, 18, 27, 
36, and 45a. The la level is 6.45 mJy beam~^. 
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Fig. 2. — Visibility amplitude versus uv distance plot for the continuum emission with 1 
a error bars. The dotted histogram is the expected amplitude for zero signal. The profile 
can be fitted with two gaussian components. The dashed line is the curve for the extended 
component and the dash-dotted line is for the compact component. The solid curve is the 
total amplitude of the two components. 
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Fig. 3. — Maps of the SiO J=8-7 emission averaged over 10 km s~^ wide intervals. The 
velocity ranges with respect to the systemic velocity of Vlsr = 5.0 km s~^ are shown in the 
upper sides of the panels. The cross marks the position of L1448C(N). The contour interval 
is 1.5 Jy beam~^ km s^^ (3 a) with the lowest contour at 1.5 Jy beam~^ km s~^ (3 a). 
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Fig. 4. — Maps of the CO J=3-2 emission averaged over 10 km s~^ wide intervals. The 
velocity ranges with respect to the systemic velocity of Vlsr = 5.0 km s~^ are shown in the 
upper sides of the panels. The cross marks the position of L1448C(N). a) The contours are 
drawn every 2.0 Jy beam~^ km s~^ (4 a) with the lowest contours at 2.0 Jy beam~^ km 
s~^ (4 a). b)-h) The contours are drawn every 1.5 Jy beam~^ km s"^ (3 a) with the lowest 
contours at 1.5 Jy beam~^ km s~^ (3 a). 
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Fig. 5. — High-velocity SiO {left) and CO (right) emission. The velocity ranges are ±51-70 
km s~^ with respect to the systemic velocity. The contours are drawn every 2.01 Jy beam~^ 
km s^^ (3cr) with the lowest contour level at 2.01 Jy beam~^ km s~^ (3 a). Green cross in 
each panel denotes the position of the continuum peak for L1448C(N). Black crosses mark 
the positions of the knots. Dashed straight line and dash-dotted line indicate the axes of BI 
(P.A. = —25°) and RI (P.A. = —20°) component, respectively. Thin solid lines are the axes 
of BII (P.A. = +15°) and RII (P.A. = -5°) components. 
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Fig. 6. — (a) The CO J=3-2 jet in the velocity range from —70 to —61 km s~^ with respect 
to the cloud systemic velocity. The cross marks the position of L1448C(N). The map is 
rotated by 25° counterclockwise. The contours are drawn every 1.5 Jy beam^^ (3 a) with 
the lowest contour of 1.5 Jy beam^^ (3 a) . (b) The CO J=3-2 jet in the velocity range 
from +61 to +70 km s~^. The map is rotated by 20° counterclockwise. The contour levels 
are the same as (a). 
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Fig. 7. — Position-velocity (P-V) diagrams of the SiO {left) and CO {right) along the jet 
axes (—25° in the blueshifted part and —20° in the redshifted part). Contours are drawn 
every 3.0 K with the lowest contour at 3.0 K in the SiO, and every 3.5 K with the lowest 
contour at 3.5 K in the CO. The horizontal thin dashed lines and thin dash-dotted lines are 
the positions of L1448C(N) and knots, respectively, while the vertical thin dashed lines label 
the systemic velocity, V^ys =5.0 km s^^. The thick line in the upper right panel denotes the 
model curve produced by C =0.6 arcsec"^ and Vq =5.0 km s^^ arcsec^^. The thick line in 
the lower right panel denote the model curve produced by C =0.8 arcsec^^ and fo =5.0 km 
s~^ arcsec"^. 
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Fig. 8. — Outflow shell model curves overlaid onto the maps of the CO J=3-2 emission 
integrated over the velocity ranges of ±1-40 km s~^. Contours are drawn every 3.8 Jy 
beam~^ km s~^ (4cr) with the lowest contour of 3.8 Jy beam~^ km s~^. a) Blueshifted 
component with a model curve (green line) with C =0.6 arcsec"^ and Vq =5.0 km s^^ 
arcsec^^. b) Redshifted component with a model curve (green line) with C =0.8 arcsec"^ 
and Vq =5.0 km s^^ arcsec"^. Both curves were produced by assumin g the inclination angle 
of the axis from the plane of the sky to be 21° (JGirart fc Acordll200ll ). 
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Fig. 9. — Outflow shell model curves in the different velocities overlaid onto the maps of 
the blueshifted CO J=3-2 at 10 km s~^ wide intervals. The maps are rotated by 25° 
counterclockwise. The model curve of each velocity was calculated using the same parameters 
as those of Fig. 8a. 
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Fig. 10. — Position- velocity (P-V) diagrams of the CO (black contours) and SiO (color 
contours) along the lines perpendicular to the jet axes. The cuts are across the emission 
peaks of knots (a) Bla, (b) Bib, c) RIa, and (d) Rib. Contours are drawn every 3.0 K with 
the lowest contour at 3.0 K in the SiO, and every 3.5 K with the lowest contour at 3.5 K in 
the CO. The vertical thin dashed lines label the systemic velocity, Kjys =5.0 km s~^. 
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Fig. 11. — Map of the CO J=3-2 emission centered at L1448C(S). The velocity offsets 
from the cloud systemic velocity are from —13 to —2 km s~^ for the blueshiftd component, 
and from +2 to +7 km s"^ for the redshifted component. Two crosses mark the positions 
of ptotostars; the southeastern one is L1448C(S) and the northwestern one is L1448C(N). 
Contours are drawn every 3 a (1.56 Jy beam~^ km s~^ in the blue and 1.11 Jy beam~^ km 
s~^ in the red) intervals with the lowest contours at 3 cr. Contours for higher flux in the 
L1448C(N) outflow are eliminated in order to provide the better view of the fainter NE-SW 
outflow. 
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Fig. 12. — SED of the continuum emission from the compact component of L1448C(N). 
The filled circles are the obser ved flux values at 3 50 GHz (t his work) , 345 GHz 
J.Tc^rgensen et ahlbnOTh. 23 GHz Jjorgensen et al.lboOTh. 1 15 GHz JSachiller et all Il995l) 
^_87 GHz JGuilloteau et aDll992h . 15 GHz JGuriel et al.lll99o] ). and 8.3 GHz fc.eipurth et al. 
2OO2I ) . The solid line shows a fit to the observed flux data using Fcxi^", where a = 1.98 
is the spectral in dex. The do t ted li ne is the expected free-free emission calculated us- 
ing the model of ICuriel et al.l (119901 ). The dash-dotted curve is the optically thick fit 
{S^ = n^B^l - exp(-r^)]) for T = 40 K. 
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Fig. 13. — Variation of the velocity centroid of the SiO emission in the redshifted part of the 
jet along the axis (P.A. -20°). 



